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ABSTRACT 

The question of the transition to global isotropy from our anisotropic local universe 
is studied using the Union 2 catalogue of Type la supernovae (SNela). We construct 
a 'residual' statistic sensitive to systematic shifts in their brightness in different di- 
rections and use this to search in different redshift slices for a preferred direction on 
the sky in which the SNela are brighter or fainter relative to the standard ACDM 
cosmology. At low redshift (z < 0.05) we find that an isotropic model such as ACDM 
is barely consistent with the SNela data at 2-3 a. A maximum likelihood analysis of 
peculiar velocities confirms this finding — there is a bulk flow of 260 kms -1 extending 
out to z ~ 0.06, which disagrees with ACDM at 1-2 a. Since the Shapley concentra- 
tion is believed to be largely responsible for this bulk flow, we make a detailed study 
of the infall region: the SNe la falling away from the Local Group towards Shapley are 
indeed significantly dimmer than those falling towards us on to Shapley. Convergence 
to the CMB rest frame must occur well beyond Shapley (z > 0.06) so this low redshift 
bulk flow will systematically bias any reconstruction of the expansion history of the 
universe. At higher redshifts z > 0.15 the agreement between the SNela data and 
the ACDM model does improve, however, the sparseness and low quality of the data 
means that the latter cannot be singled out as the preferred cosmological model. 

Key words: cosmology: theory, dark matter, large-scale structure, peculiar velocities, 
cosmic microwave background, cosmological parameters, Type la supernovae 
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Modern cosmology is founded on the cosmological princi- 
ple which assumes that the universe is homogeneous and 
isotropic. The local universe is however observed to be 
anisotropic and inhomogeneous, exhibiting the 'cosmic web' 
of voids and superclusters. This presumably causes the Lo- 
cal Group of galaxies moves towards a preferred direction 
i = 276° ± 3,6 = 30° ± 2 at 627 ± 22 kms"\ as inferred 
from the dipole anisotropy of the cosmic microwave back- 
ground (CMB) radiation ( |Kogut et al.|1993 l. 

On the other hand the overall isotropy of the CMB 
(barring the dipole anisotropy due to our local motion) 
provides strong evidence for an isotropic universe on very 
large scales^] Where does the transition between these two 



1 The WMAP observations of anomalies in large-angle 
anisotropies in the CMB, e.g. the hemispherical asymmetry ( |Erik-| 
|sen et aL"||2004[l and the unexpected quadrupole-octupole align- 
ment ( |de Oliveira-Costa et al.|2004[ |, have led many to question 
whether the CMB is indeed statistically isotropic (e.g. |Copi et al.| 



regimes occur? While high-quality data exist at low redshifts 
and the CMB provides reliable information at very high red- 
shifts, the data is rather sparse and mainly of poor quality on 
the intermediate scales of interest. Here the SNela Hubble 
diagram is the key source of information and so we use the 



comprehensive Union 2 catalogue (Amanullah et al. 20101 



to investigate these important questions. 

SNela data has been examined previously to test the 



isotropy of the universe ( Kolatt & Lahav 



2001 



Bonvin, Dur- 



rer fe Kunz|[2006| |Gordon, Land fc Slosar 2007| |Schwarz & 

Laskarl 12008 



Weinhorst 200 7| |Gupta, Saini & Las kar||2008| |Koivisto _ 
Mota|2008a|b||Blomqvist, Enander fe Mortsell|2008||Coor ~ 
Holz fc Caldwell||2008| |Gupta fc Saini||2010| |Koivisto et al 



20101 and to determine whether the universe accelerates dif- 



ferently in different directions. Recently a marginal (la) de- 
tection of anisotropy has been reported on spatial scales 



l |2007[ l). However others have argued that these anomalies may 
simply be due to the manner in which the galactic foreground 
was masked l |Pontzen fc P ciris 2010;). We look to forthcoming 
observations by Planck to resolve this contentious issue. 
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where dark energy becomes important (Cooke & Lynden 



Bell|2010|[Antoniou fc Perivolaropoulos|2010[ >. Clearly, bet- 
ter quality and more complete surveys are needed before any 
firm conclusions can be drawn. However, although these de- 
tections are not significant by themselves, a puzzling and 
perhaps accidental feature is the alignment of the detected 
anisotropy with the the CMB dipole direction. In this work, 
we demonstrate that the alignment at low redshift is due 
to the attraction of huge structures such as the Shapley su- 
percluster. At high redshift, the alignment seems to become 
statistically insignificant but given the sparse and poor qual- 
ity data, no strong conclusions can be drawn. 

On small scales, the CMB dipole and the bulk flow are 
aligned, which is unsurprising as the common source of both 
these motions is very likely the anisotropic distribution of 
matter in the local universe. However, a bulk flow much 
larger than ex pected has been found extending out to at 
least 120 Mpc fWatkins et"aT||2009| |Lavaux et al.||2010| ), j 2 ] 
which is a ~ 2 — 3a fluctuation in a ACDM model since 
convergence to the CMB rest frame ought to occur at much 
smaller scales in this model. At low redshifts, we study the 
bulk flow using two different methods: first by a method of 
'smoothing and residuals' (Sections [3] and |4| and secondly 
by a maximum likelihood analysis ( Section (5T. We show that 
these two methods are in good agreement at small redshifts 
(z < 0.05) and confirm that there is indeed a discrepancy 
between the ACDM model prediction for the bulk flow and 
the SNela observations. 

The Shapley concentration at z ~ 0.035 — 0.055 is 
believed to be the main source of our large bulk motion. 
We study the infall region around Shapley (Section [6| and 
demonstrate that SNe la beyond Shapley are systematically 
brighter than expected in an isotropic universe (as they are 
falling towards Shapley), while SNela between the Local 
Group and Shapley are systematically dimmer (as they are 
also falling towards Shapley, but away from us). This result 
is obtained using our smoothing and residuals scheme. 

At high redshift, the \ 2 statistic cannot be used since 
the bulk flow becomes small relative to the Hubble expan- 
sion rate. Using the method of smoothing and residuals for 
z > 0.05 we find that an isotropic model such as ACDM is 
consistent with the data. However, the poor quality of the 
data means that anisotropic models cannot be eliminated. 

Since the high redshift results may be contaminated by 
the low redshift data, we perform a 'cosmic tomography' — 
the data is sliced up in redshift and the question of isotropy 
is studied separately for each slice. The differential analysis 
is then complemented by a commulative analysis to estab- 
lish the correlation between different shells. We present our 
conclusions in Section 



2 THE UNION 2 SUPERNOVAE CATALOGUE 



20091 and SDSS (Kessler et al. 20091, in addition to data 



The Union 2 catalogue (Amarmllah et al. 20101 contains 
557 SNela of which 165 are at z < 0.1. This is the largest 
compilation of SNe la to date from several different surveys 



including Union ( Kowalski et al. 2008 1 , Cf A ( Hicken et al 



2 |Kashlinsky et aL ] | |2008| |2010^ have reported an even faster 
coherent flow out to at least ~ 800 Mpc, using measurements of 
the kinematic Sunyaev-Zeldovich effect in galaxy clusters. 



on 6 new SNela reported by Amanullah et al. (20101. While 



efforts have been made to obtain good control over the sys- 
tematics, it is recognised that combining data from different 
surveys can introduce additional biases. Nevertheless this is 
the best data set available at present. 

In Figure[T] we illustrate by the size of the spots the dis- 
crepancy between the luminosity distance of Union 2 SNe la 
and the value predicted by the standard ACDM model, for 
z < 0.06 (top panel) and z > 0.06 (middle panel). The red 
spots indicate SNe la further away than the ACDM predic- 
tion, while the green spots indicate those which are closer. 
In the bottom panel of Figure [l] the size of the spots repre- 
sents the observational uncertainties. We note that there is 
a clear correlation between the latter and the discrepancy 
between the luminosity distances and the ACDM model. 

As seen in Figure [2] the data is quite homogeneously 
distributed over the sky at redshifts z < 0.06; however at 
higher redshifts the data becomes increasingly sparse. 



3 THE METHOD OF RESIDUALS 

If the assumption of the isotropy of the universe is correct 
then any residuals remaining from the subtraction of the 
isotropic model from the data should be distributed ran- 
domly around zero (assuming the systematics are under con- 
trol). In this work, we fit the standard ACDM model to the 
Union 2 data (Amanulla h et al.|2010 |, subtract the best-fit 
model from the data and then build a statistic to analyse 
the residuals. This analysis involves the following steps: 

First, we consider a H(z) parameterization (viz. 
ACDM) and a set of SNela data Hi(zi,9i,((>i) — we choose 
to use galactic coordinates (I, b) for 9i,cj>i. At low redshifts 
z < 1, the Hubble law implies a linear relationship between 
distance and redshift so the choice of cosmological model is 
irrelevant, however this becomes important at high redshift. 

Second, we obtain the best fit form of H(z) and the 
corresponding best-fit distance moduli p,(z), then construct 
the error-normalised difference of the data from the model, 



Qi(Zi, 0i, <t>i) = 



m(zi,6i,(j}i) — fjj(zi, 6j, fa) 

cn[zi, 0i, 4>i) 



(1) 



as in previous work ( Perivolaropoulos & Shafieloo 2009 
|Shafieloo, Clifton &: Ferreira||2010| >. Here, Hi(zi,9i,4>i) and 
o~i(zi, 9i,<j>i) represent the distance modulus and the relative 
error bar of the i'th data point, and jii(zi, 8i, <j>i) is the ex- 
pected distance modulus in the assumed theoretical model 
at Zi. Henceforth we work with these residuals, qi(zi, 6i,4>i), 
and consider their angular distribution on the sky. 

Third, we define a measure Q(6, <f>) on the surface of a 
sphere of unit radius using these residuals: 



(2) 



where N is the number of SNela data points and 
W(9, <j>, 9i, <pi) is a weight function (or window function) in 
our 2D smoothing over the surface of the sphere. We define 
this weight by the Gaussian distribution: 



i 



2nS 



exp 



Lie, 



2<5 2 



(3) 
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Figure 1. Top panel: Aitoff-Hammer presentation of the Union 2 
data for z < 0.06. The red spots represent the data points with 
distance moduli, /^data; larger than the values, ^tACDMi predicted 
by ACDM and the green spots are those with /^data < A*ACDm; 
the spot size is a relative measure of the discrepancy. A dipole 
anisotropy is visible around the direction b = — 30°, = 96° (red 
points) and its opposite direction b = 30° , i = 276° (small green 
points), which is the direction of the CMB dipole. The middle 
panel is the same plot for z > 0.06. The data seems to be homo- 
geneously distributed at small redshifts but suffers from a signif- 
icant selection effect at large redshifts. The bottom panel shows 
the full Union 2 data set (0.015 < z < 1.5) but now the size of 
the spots correspond to the observational uncertainties — a clear 
correlation is seen with the model-data discrepancy. 
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Figure 2. The distribution in galactic longitude (£) and latitude 
(ft) of the Union 2 SNela, as a function of the redshift z. At 
high redshift, there are significant selection effects and it is diffi- 
cult to draw any definite conclusions about the anisotropy of the 
universe; however at small redshift (z < 0.06) the data is quite 
homogeneous and robust tests of local anisotropy can be made. 



where 8 is the width of smoothing and L(9,<f>,9i,(f>i) is the 
distance on the surface of a sphere of unit radius between 
two points with spherical coordinates (6, cf>) and (9i,(f>i): 

R 



L(6, 4>, 0i, 4>i) = 2 arcsin 



2 ' 



(4) 



where 



R 



>in(#i) cos(0i) — sin(#) cos(<j!>)] + 

2 2\ 1//2 
[sin(#i) sin(i^i) — sin(#) sin(</>)] + [cos(#i) — 008(6*)] 

Thus, any anisotropy in the data in any specific direction 
will translate into Q(8, <f>) being significantly less or more 
than zero, depending on the quality of the data. 

Finally we adopt a value for S, calculate Q(8, <j>) on the 
whole surface of the sphere and find the minimum and maxi- 
mum of this quantity, Q(6> min , min ) and Q(6 max , max ). Our 
statistical measure of anisotropy is based on the difference 



Q(6 n 



(5) 



i.e. a large value of AQdata implies significant anisotropy. 
To estimate the significance, we simulate 1000 realizations 
of the data with the same angular positions on the sky and 
the same error bars (assuming a Gaussian distribution), de- 
termine A Qj and then do a simple test to determine the 
significance of our results for the real data. Any anisotropy 
in the SNe la sample would be more significant if A Qdata is 
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larger than most of A QjS for the simulated data. From the 
resulting frequency distribution of this statistic one can de- 
rive a P-value, defined as the probability that, given the null 
hypothesis, the value of the statistic is larger than the one 
observed. We remark that in defining this statistic one has 
to be cautious about a posteriori interpretations of the data; 
a particular feature observed in the real data may be very 
unlikely, but the probability of observing some feature may 
|Spergel et al.| fl2003|), |Lewis] (|2003|) a nd 



be quite large — see 



the discussion in Hamann, Shafieloo & Souradeep (20101 



We would like to point out here some additional ad- 
vantages of our method. First, the only parameter in this 
analysis is the value of 5 — the width of gaussian smooth- 
ing on the surface of the 2D sphere. Its value cannot affect 
the results when testing the consistency of the data with an 
isotropic universe, since we are dealing only with the residu- 
als and the real data is compared with its simulated samples 
in exactly the same manner. Second, the isotropy of the data 
can be checked in a completely model-independent manner 
and the significance of the findings can also be checked to 
avoid any misinterpretation of the data. Third, the method 
is able to detect anisotropy even in small patches of the sky 
if there is sufficient data available. In fact our method is not 
limited to finding only dipoles in the sky but can be useful if 
there is a neighboring local void or any unexpected bulk flow 
up to modest redshifts in a specific direction. Moreover, be- 
cause of working with the normalised residuals, the method 
is insensitive to those regions of the sky where there is no 
data, which makes it relatively bias-free. Hence it is partic- 
ularly suitable for z > 0.05 where there are rather few data 
points. Although in this analysis we focus on the isotropic 
model of the universe, our method can be used to test any 
other model, in particular anisotropic models. 

We wish to emphasise the difference between testing the 
consistency of a cosmological model with the data, and find- 
ing the model that best describes the universe. Depending 
on the quality, quantity and distribution of the data there 
may be several (degenerate) models that are all consistent 
with the data. For example, we may find that the isotropic 
model of universe fits the data, but that the data is so sparse 
at high redshifts that in many patches of the sky we do not 
have even a single data point, so an anisotropic model which 
has features in these directions can also fit the data. 



4 TEST OF ISOTROPY USING THE METHOD 
OF RESIDUALS 

We now use the method of smoothing and residuals to test 
the isotropy of the universe using SNela data at different 
redshifts and for different values of the smoothing parame- 
ter 8. We divide the Union 2 data ( Amanullah et al.]|2010[ ) 
into 5 bins at z < 0.1 and 7 bins at z > 0.1 so as to have 
sufficient number of data points in each bin — see Table [l] 
In each case the real data is compared with the results from 
1000 Monte Carlo realisations assuming an isotropic uni- 
verse. The underlying cosmological model is assumed to be 
a flat ACDM model with f2om = 0.27 (which is very close to 
the best-fit flat ACDM model for Union 2 data). 

In Figure [3] we show results for the first 3 low redshift 
shells and the cumulative result when data in all 3 shells 
is combined. Although each individual shell is reasonably 



consistent with isotropy, combining the data indicates an 
inconsistency of the ACDM model with the data at > 2a. 
This happens because the direction of the mild anisotropy 
in all 3 shells point the same way. 

We have added Table [2] to emphasise that the impor- 
tant quantity at low-z is not the P-value alone, since similar 
P-values are obtained also at high redshift. One should con- 
sider also the consistency of adjacent shells: the anisotropy 
in the first 3 shells point towards nearly the same direction. 
This strong correlation makes the P-value for the cumulative 
result of the first 3 shells rather low: z = 0.045. Figure [3] and 
Table [2] show that at z > 0.1 there are no adjacent shells 
that point towards the same direction, hence, the effect of 
these shells at high redshift cancel each other in the cumu- 
lative result and the effect of first 3 shells dominates. The 
importance of our differential shell analysis is now evident 
— any anisotropy in some specific redshift range would be 
smoothed out if one were to use only the whole dataset. 

It is also possible to look for dipoles in different shells 
using our method. We look for a direction in the sky for 
which the value of AQ between that point (I, b) and its op- 
posite point on the surface of the sphere — i.e. (£+ 180, —6) 
if £ < 180 and (£ - 180, -6) if I > 180 — is maximum. 
To calculate dipoles in different shells we set S = n/2; for 
this value the weighted average using gaussian smoothing 
is almost identical with the weighted average by the inner 
product (i.e. cos 9) that is usually used to find dipoles. 

In Figure [4] we see the positions of the dipoles in galac- 
tic coordinates for different redshift shells. The size of the 
spots are proportional to the P-values that are derived by 
comparing the results from the real data to the 1000 Monte 
Carlo realizations of the data when the assumed model is 
ACDM. Bigger spots indicate a more anisotropic behavior 
or lower level of agreement between ACDM and the real data 
(the size of the spots is proportional to the P-value). In the 
bottom panel we see the cumulative results which can help 
us to identify correlations between neighboring shells. 

In Table [2] we see these results for different shells, both 
individual and cumulative. It is interesting to see that the 
P-value drops to 0.014 at z < 0.045 i.e. a clear inconsistency 
between ACDM (i.e. isotropy) and the data. The P-value 
increases again with increasing redshift, suggesting that we 
are approaching the CMB rest frame. This can be seen more 
clearly by looking at the top panel of Figure[5]where P-values 
are plotted versus redshift for different cumulative shells. In 
the bottom panel of Figure [5] we can see the results from 
our residuals analysis where we have divided the data into 
two subsamples — one for low redshifts 0.015 < z < 0.06 
and another for higher redshifts 0.15 < z. At small redshifts, 
the isotropic universe lies over a away from the data with 
P=0.054 for 0.015 < z < 0.06. However at higher redshift 
this discrepancy drops to ~ la with P=0.594 for z > 0.15. 



5 TEST OF ISOTROPY AND LOCAL BULK 
FLOW USING MAXIMUM LIKELIHOOD 
ANALYSIS 

Now we use a different method to examine the low-redshift 
aniosotropy detected above. We assume a bulk flow and use 
maximum likelihood analysis method to find its value. We 
consider SNela at z < 0.06 (since the bulk flow would be 
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Table 1. Number of SNela per redshift bin in the Union 2 catalogue 



Az 


0.015-0.025 


0.025-0.035 


0.035-0.045 


0.045-0.06 


0.06-0.1 


0.1-0.2 


N 


61 


48 


18 


15 


23 


55 


Az 


0.2-0.3 


0.3-0.4 


0.4-0.5 


0.5-0.6 


0.6-0.8 


0.8-1.4 


N 


62 


63 


58 


43 


51 


60 



< 



Union 2 compilation 
68% and 95% CL 



0.015 < z < 0.025 



< 



Union 2 compilation ■ 
68% and 95% CL ■ 




0.025 < z < 0.035 




Figure 3. The anisotropy measure AQ (eq|5jl plotted against the smoothing length S, for different redshift shells. The full (red) lines 
show the results for the Union 2 SNela data and the dotted (blue) lines are 68% and 95% confidence limits derived from 1000 Monte 
Carlo realisations. While each individual shell is consistent with isotropy, all shells have a mild anisotropy in the same direction so the 
cumulative result from combining them is inconsistent with isotropy at over 2a, as seen in the lower right panel. 



negligible in relation to the Hubble expansion rate at higher 
redshift) and minimise, for the Union 2 data: 

X 2 (Qom,V bulk ,j(£, b)) = 

N r,,Union2 K1 /jACDM m _ , V hulk .j(e,b) \ ggj 2 



E 



51og 10 (d^ UM (n 0m ,H ) + 



where N is the number of SNela, j(£,b) is a unit vec- 
tor representing the direction on the sky in Galactic co- 
ordinates, Vbuik is the local bulk flow velocity in kms -1 , 



d£ CT>M (Qom, Ho) is the luminosity distance for the standard 
flat ACDM model and Qom is the present day matter den- 
sity (set equal to 0.27 although the choice is not important 
since the Hubble law is linear here). It should be noted that 
(ffi dealing with SNela data, Hq is treated as a 'nuisance 
parameter'. We consider different cumulative redshift slices 
of the data, viz. 0.015 < z < 0.025, 0.015 < z < 0.035, 
0.015 < z < 0.045 and 0.015 < z < 0.06. 

In Figure [6] we see that the data for 0.015 < z < 0.025 
(top panel) fits best with Vbuik = 260 kms -1 in the direc- 
tion b = 16°, I = 271° which is close to the direction of the 
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Table 2. Results from the residuals analysis for the P-value which quantifies the level of agreement between the isotropic universe and 
the data, for shells in redshift. On the left are the results for the individual shells and on the right for the cumulative shells. 



Az 


b° 


e° 


P-value 


Az 


b° 


t° 


P-value 


0.015-0.025 


46 


265 


0.140 


0.015-0.025 


46 


265 


0.140 


0.025-0.035 


2 


320 


0.112 


0.015-0.035 


27 


297 


0.039 


0.035-0.045 


20 


312 


0.354 


0.015-0.045 


26 


300 


0.014 


0.045-0.06 


-46 


65 


0.267 


0.015-0.06 


19 


309 


0.054 


0.06-0.1 


-41 


75 


0.637 


0.015-0.1 


14 


316 


0.153 


0.1-0.2 


-3 


219 


0.608 


0.015-0.2 


13 


299 


0.265 


0.2-0.3 


29 


21 


0.602 


0.015-0.3 


20 


314 


0.226 


0.3-0.4 


-83 


348 


0.854 


0.015-0.4 


10 


314 


0.329 


0.4-0.5 


70 


238 


0.428 


0.015-0.5 


27 


307 


0.285 


0.5-0.6 


16 


15 


0.177 


0.015-0.6 


27 


326 


0.180 


0.6-0.8 


-77 


45 


0.108 


0.015-0.8 


3 


332 


0.279 


0.8-1.4 


-54 


152 


0.947 


0.015-1.4 


-2 


332 


0.369 



0.4< i. <0.5 
O 





CMB dipole (b = 30 V = 276°) fKogut et aL]p93] ). How- 
ever, the la contours are quite large and an isotropic local 
universe with Vbulk = is also consistent with the data at 
2a. When we assume a bigger cumulative cut, the la and 2a 
regions shrink and we can see that the isotropic universe is 
now shifted outside the 2a region (Figure [6] bottom panel). 
We can observe the same trend up to 0.015 < z < 0.045 
where the 1 and 2 a regions are still shrinking as seen in 
Figure [7] (top panel). However in the next step when we 
consider the data for 0.015 < z < 0.06, the 1 and 2 a re- 
gions start becoming larger again (Figure [TJ bottom panel). 
These results are completely consistent with our results from 
the smoothing and residuals analysis. So we can conclude 
that the bulk of the local universe of radius z ~ 0.05 is 
moving towards the CMB dipole direc tion with Vbulk = 260 
kms" . Our results are consistent with 
who estimate a bulk flow of Vbulk 
wards b = 60° (+6, -6),£ = 282°(+ll 



Watkins et al. 



(20091 

416 ± 78 kms^to- 
-11) at a scale of 
100ft" 1 Mpc. On this scale we find Vbulk = 250 (+190, -160) 
kms" 1 towards b = 21°(+34, -52), £ = 287° (+62, -48), i.e. 
our results are consistent within la. 

In Figure [8j we see that the derived bulk flow using the 
cumulative SNe la data is inconsistent with the prediction of 
a flat ACDM model at 1-2 a. This is a slight improvement 
over using residuals (where the disagreement with ACDM 
was at 2-3 a) the reason being that at small redshift one 
needs to compare the data with the model at first-order in 
perturbation theory. At high redshifts the comparison can be 
done directly between the data and the unperturbed model 
since the perturbations have subsided on large scales. 



Figure 4. Top panel: The dipole direction found by the residuals 
method in different redshift shells. Bottom panelThc cumulative 
dipole direction in shells of increasing radii. The size of the spots 
is proportional to the P-value; larger spots represent more signif- 
icant disagreement between the data and the ACDM model. 



6 SHAPLEY INFALL 

It is generally believed that our local bulk flow is due mainly 
to the gravitational attraction of the Shapley supercluster 
which lies at around z ~ 0.04 and that beyond Shapley the 
reference frame should converge to the CMB rest frame. We 
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Figure 5. Top panel: P-value for the consistency of the isotropic 
universe with the data versus redshift. At z 0.05 the P-value 
drops to 0.014 showing that isotropy is excluded at 3<r and we are 
not yet in the CMB rest frame. Bottom panel: The cumulative 
analysis shows that at small redshift isotropy is excluded at 2- 
3 a with P = 0.054 for 0.015 < z < 0.06; however at higher 
redshift agreement is achieved within 1 a, with P = 0.594 for 
0.15 < z < 1.4. 



use the method of smoothing and residuals, discussed in Sec- 
tion]^] to study the infall towards the Shapley concentration. 
Somewhat surprisingly the SNela data provides us with a 
clear picture of the infall. 

We take Shapley to have an approximate extension of 
0.035 < z < 0.055. First we consider all data in the redshift 
band 0.015 < z < 0.0345 which contains 109 SNela and 
evaluate the dipole which, as shown in the Figure [9j points 
towards Shapley (which is approximately aligned with the 
direction of the CMB dipole) . While the movement towards 
Shapley seems to be strongly favored by the data, the P- 
value for the isotropic universe is about P — 0.039. Next, we 
consider the data in the redshift range 0.0522 < z < 0.095 
which contains 32 SNela and as shown in Figure [9] the di- 
rection now becomes completely opposite, indicating these 
SNela are infalling towards Shapley. The P-value for the 
isotropic universe in this range is P = 0.339. Thus future 
precision SNe la data can provide us with valuable informa- 
tion on the formation history of this giant supercluster. 




Figure 6. The top panel shows the dipole from the maximum 
likelihood analysis in the redshift band 0.015 < z < 0.025 (61 
SNela). The best fit point is at (b = 16°, £ = 271°) for V hulk = 
250 kms -1 and the red and green contours are the 1 and 2 a 
confidence regions. The large blue spot is the direction of CMB 
dipole (b = 30° ± 2,1 = 276° ± 3). The larger yellow spot, close 
to the CMB direction, is the best fit direction from the residuals 
analysis and the smaller black spot is the best fit direction from 
the maximum likelihood analysis (the magnitude of the dipole 
was given in Figure [sj . The bottom panel shows the dipole for 
the redshift range 0.015 < z < 0.035 (109 SNela). The blue spot 
is the CMB dipole and the black spot at b = 21°, £ = 287° is the 
best fit from the likelihood analysis for Vbulk = 260 kms -1 , while 
the red and green patches are the 1 and 2 a confidence regions. 

7 CONCLUSIONS 

We have used SNela data from the Union 2 catalogue to 
study the (an)isotropy of our universe. Since the low and 
high redshift anisotropy can in principle be of different ori- 
gin, we have analysed the data tomographically in different 
redshift slices, and also the cumulative data set in order to 
trace any correlation between different slices. We have devel- 
oped a statistical tool of 'residuals' to examine the data in 
an unbiased manner at all redshifts. At low redshift, we have 
performed a maximum likelihood analysis (which is however 
unsuitable at high redshift). We find that an isotropic model 
like ACDM is 2-3 a away from the data at z < 0.06. Al- 
though the agreement improves at high redshifts (to within 
1 a), we cannot single out ACDM as the preferred model of 
the universe. The data becomes rather sparse at high red- 
shift and the error in distance measures increases, so the 
data may also agree with alternative anisotropic models. 
At low redshift, our results are rather robust and we find 
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Figure 8. The bulk flow as a function of redshift from the like- 
lihood analysis. We see that a fast flow with Vb u ik = 260 kms — 1 
persists up to at least z = 0.06 and systematically exceeds the pe- 
culiar velocity (y3 times the 1-d RMS value) expected in ACDM 
normalised to WMAP-5 parameters (Wat kTns et al.| |2009). Note 
that the 1-d RMS expectation is shown (blue line) in order to 
emphasise that statistical analyses are carried out on the sepa- 
rate components of the velocity, and also to facilitate comparison 
with the previous result of |Watkins et ah] | |2009| . 




Figure 7. The top panel is for the range 0.015 < z < 0.045 (127 
SNela) and the bottom panel is for 0.015 < z < 0.06 (142 SNela). 
The blue spot is the CMB dipole (b = 30° ± 2, £ = 276° ± 3) and 
the black spots are the best-fit from the likelihood analysis at 
(b = 15°, I = 291°) for Vb ulk = 270 kms" 1 for the top panel 
and at (6 = 8°,£ = 298°) for V bulk = 260 kms" 1 for the bottom 
panel, while the red and green patches are the 1 and 2 a confidence 
regions. 



We also note the interesting observation by |Tsagas| 



(2010) that observers with peculiar velocities have local 



expansion rates which are appreciably different from the 
smooth Hubble flow, so can experience apparently acceler- 
ated expansion when the universe is actually decelerating. 
Thus whether dark energy really needs to be invoked to 
explain the SNela Hubble diagram and other observations 
remains an open question ( Sarkar| 2008). 



a bulk flow of about 260 Km s _1 in the direction of the Shap- 
ley supercluster. We show that the Union 2 data provide the 
first evidence of the infall on to Shapley; SNefa which are 
falling away from us and towards Shapley are statistically 
dimmer that those which lie beyond this supercluster and 
are falling towards us. We see no indication of the decay of 
the bulk flow after Shapley which suggests that the scale 
of anisotropy of our local universe is bigger than is usually 
assumed and extends beyond z ~ 0.1. 

Our analysis and results are important for the study 
of the expansion history of the universe and the properties 
of dark energy. In all SNela compilations, an uncertainty 
of 300-500 Kms" 1 is assumed for each data point to allow 
for bias introduced by random peculiar velocities. However 
when there is a coherent motion of SNe la towards a specific 
direction, this bias cannot be removed by just increasing 
the size of the error bar (i.e. assuming the peculiar veloci- 
ties to be random). We will present in future work the ef- 
fect of this systematic motion of SNe la at low redshifts on 
the reconstruction of the expansion history of the universe 
and estimation of cosmological parameters like q(z), w(z) 
or 'Om'(z) (|Sahni, Shafieloo fc Starobinsky|2008[|Shafieroo,| 
Sahni fc Starobinsky|2009[ ). 
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